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Abstract .  This  paper  is  concerned  with  conduction  of  heat  and  the  related 
problem  of  propagation  of  thermal  waves  in  stationary  rigid  solids . 

Special  attention  is  given  to  rate-dependent  response  and  the  ensuing 
conditions  of  propagation  in  the  conducting  medium.  By  considering  small 
time -dependent  temperature  variations  superposed  on  a  finite  nonuniform 
equilibrium  temperature  field,  certain  conclusions  are  reached  which  also 
shed  light  on  the  so-called  "second  sound"  phenomenon. 


1 .  Introduction 

This  paper  is  concerned  with  heat  conduction  and  propagation  of 
thermal  waves  in  rigid  solids.  Our  approach  and  point  of  view  is  that  of 
recent  developments  in  continuum  thermodynamics.  We  employ  the  balance 
of  energy  and  the  (Clausius-Duhem)  entropy  production  inequality,  as  well 
as  constitutive  assumptions  characterizing  the  local  thermal  behavior  of 
the  material.  For  example,  if  temperature,  temperature  gradient  and  rate 
of  temperature  are  taken  as  independent  variables,  then  the  thermal 
variables  such  as  specific  internal  energy,  specific  entropy  and  heat 
flux  are  determined  by  constitutive  equations.  When  the  constitutive 
assumptions  include  rate  of  temperature  as  an  independent  variable  there 
arises  a  local  production  of  entropy  (see  Eq.  (2.2))  in  addition  to  the 
entropy  production  due  to  conduction,  which  is  inherent  in  the  usual  rate- 
independent  response.  The  rate-dependent  thermal  response  studied  here 
and  the  corresponding  generalization  of  the  classical  (Fourier)  heat 


conduction  equation  leads  to  a  prediction  of  finite  speed  thermal  waves, 
as  well  as  to  other  novel  phenomena  associated  with  local  thermal 
dissipation . 

Before  describing  the  scope  of  the  paper  in  detail,  it  is  desirable 
to  recall  for  background  information  certain  other  developments  which 
bear  on  conduction  of  heat  in  rigid  solids.  Within  the  framework  of 
modern  developments  in  continuum  physics,  attention  was  largely  con¬ 
fined  until  recently  to  purely  mechanical  constitutive  equations.  Argu¬ 
ments  for  obtaining  necessary  and  sufficient  conditions  for  the  validity 
of  the  entropy  production  inequality,  limited  to  nonlinear  elastic 
materials  with  heat  conduction  and  viscosity,  were  given  by  Coleman  and 
Noll  C l] .  A  more  extensive  discussion  of  generalizations  of  the  classical 
theory  of  linearly  viscous  fluids  with  linear  heat  conduction,  in  which 
the  independent  thermodynamic  variables  are  temperature  and  temperature 
gradient,  is  contained  in  a  paper  by  Coleman  and  Mizel  [2].  Among  subse¬ 
quent  developments,  we  mention  Coleman’s  work  [3]  on  thermodyanmies  of 
materials  with  fading  memory,  in  which  the  thermo-mechanical  constitutive 
equations  are  assumed  to  be  functionals  over  the  time  histories  of  the 
chosen  independent  variables.  Similar  constitutive  assumptions  are  uti¬ 
lized  bjr  Coleman  and  Gurtin  [4]  in  their  investigation  of  the  thermal 
behavior  of  rigid  heat  conductors. 

Certain  types  of  well-known  mechanical  behavior  -  typical  of  materials 
or  rate-type  —  such  as  that  for  linear  (Newtonian)  viscous  fluids  are 
appropriately  characterized  by  rate -dependent  functions  rather  than  by 


memory  functionals.  Likewise,  we  may  expect  certain  thermal  behavior  of 
rigid  conductors  to  be  best  described  by  rate-dependent  functions  and 
this  is  the  underlying  premise  of  the  generalization  of  the  classical 
heat  conduction  equation  sought  here. 

That  a  generalization  of  the  classical  heat  conduction  equation  for 
solids  is  desirable  has  been  recognized  for  several  years,  This  view  is 
based  primarily  on  the  fact  that  the  classical  equation  is  parabolic  and 
does  not  admit  the  possibility  of  finite  speeds  of  propagation  of  thermal 
pulses.  Most  of  the  previous  attempts  to  alter  the  heat  conduction  equa¬ 
tion  so  as  to  predict  thermal  waves,  often  referred  to  in  the  literature 
as  "second  sound,"  have  been  based  on  an  ad  hoc  generalization  of  Fourier's 
heat  conduction  law;  see  for  example  Chester  [5]  and  Ulbrich  [6] .  This 
generalization,  referred  to  by  Ulbrich  as  Vernotte's  hypothesis  [7], 
appears  as 

q  +  t  q  ■  k  grad  0  (l.l) 


We  may  recall  here  that  constitutive  equations  of  linear  viscous  fluids 
can  be  regarded  as  exact  constitutive  relations  which  describe  the 
behavior  of  a  class  of  fluids  in  all  motions.  On  the  other  hand,  in 
a  functional  theory  (e.g.,  of  the  type  developed  by  Coleman  [3^)j  a 
linearly  viscous  material  approximates  a  general  material  with  fading 
memory  only  in  the  limit  of  "slow"  motions. 
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for  rigid  conductors,  where  q,  ij,  8  represent  the  heat  flux,  the  time 
rate  of  heat  flux,  and  temperature,  while  t,  k  are  material  constants. 
With  reference  to  (l.l)  Ulbrich  [6]  states  "Although  Vernotte's  proposed 
revision  of  Fourier's  hypothesis  adequately  circumvents  the  paradox  of 
infinite  velocity,  no  apparent  physical  justification  can  be  offered  for 
the  addition  of  the  second  term  ..." 

A  different  approach  from  that  based  on  (1.1),  but  with  the  same  ob¬ 
jective,  is  taken  by  Kaliski  [83.  He  starts  with  the  assumption  that  the 
heat  conduction  equation  should  be  a  second  order  hyperbolic  equation 
and  then  generalizes  Onsager's  reciprocal  relations  and  the  associated 
form  of  the  entropy  inequality  so  as  to  accommodate  this 
assumption. 

We  also  draw  attention  to  a  very  recent  paper  by  Gurtin  and  Tipkin  [9] 
whose  work  is  motivated  by  (1.1).  They  develop  a  nonlinear  theory  for 
rigid  heat  conductors  in  which  the  constitutive  assumptions  are  in  the 
form  of  functionals  over  the  temperature  history.  They  show  that  an 
equation  of  the  form  of  (l.l)  is  a  special  case  of  their  linearized 
constitutive  equation  for  the  heat  flux  q  .  Their  associated  general- 

A* 

ized  heat  conduction  equation  in  the  linear  theory  is  given  by 


20(x,t)  +  p(0)fi(x,t)  +  I  0'(s)e(x,t-s)ds 

^  A*  V  q  ^ 

■  a(O)A0(x,t)  +  J  a'(s)&0(x,t-s)ds  +  r(x,t) 

A#  V  A  A  A 


(1.2) 


4. 


'-■'MB' 


where  a(s),  fj(s)  are  called  the  heat  flux  relaxation  function  and  the 
energy  relaxation  function.  This  equation  is  in  a  form  that  predicts 
finite  speeds  of  propagation  of  thermal  waves.  An  unusual  feature  of 
(1.2),  however,  is  the  presence  of  r  (rather  than  r)  which  evidently  is 
due  to  the  fact  that  (1.2)  is  obtained  from  the  first  time  derivative  of 
the  energy  equation  rather  than  from  the  energy  equation  itself. 

In  the  present  paper,  which  is  concerned  with  rigid  stationary  con¬ 
ductors,  in  general  we  use  direct  (coordinate  free)  notation.  Thus,  vectors 
and  points  in  three-dimensional  Euclidean  space  are  denoted  by  boldface 
Latin  lower  case  letters  while  boldface  Latin  capital  letters  are  used  to 
designate  tensors  of  order  two;  also,  Greek  lower  case  letters  are  reserved 
for  scalar  thermodynamic  variables.  An  exception  to  the  above  notation 
occurs  when  we  write  the  energy  equation  as  a  partial  differential  equa¬ 
tion  in  the  temperature.  Then,  for  ease  of  comparison  and  analysis,  all 
quantities  in  the  energy  equation  are  expressed  in  terms  of  their  Cartesian 
components . 

In  section  2,  following  seme  preliminaries,  for  clarity  and  later 
comparison  we  include  a  discussion  of  nonlinear  constitutive  equations 
for  rate-independent  response  when  temperature  and  temperature  gradient 
are  taken  as  independent  variables.  A  generalization,  in  which  the  inde¬ 
pendent  variables  include  the  rate  of  temperature,  is  introduced  and 
developed  in  section  3«  This  leads  to  a  local  thermal  dissipation  of 
energy  and  to  the  possibility  of  temperature  changes  without  heat  flow. 

The  appropriate  "heat  conduction  equation"  is  then  derived  and  is  used 

These  constitutive  equations  when  properly  linearized  reduce  to  those 

appropriate  to  the  classical  theory  of  heat  conduction. 
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in  section  4  to  investigate  the  existence  of  finite  speeds  of  propagation 
of  second  order  discontinuities  in  the  temperature. 

In  order  to  arrive  at  more  explicit  conclusions,  we  special^-e  in 
section  5  the  results  of  section  3  by  assuming  that  the  dependence  on 
the  temperature  rate  is  linear.  We  then  re-examine  the  energy  equation 
for  the  case  of  isotropic  materials  with  a  center  of  symmetry  and,  in 
particular,  make  certain  observations  regarding  the  speed  of  propagation 
of  thermal  waves  in  one  space  dimension.  Of  special  significance  is  the 
result  in  section  5  that  real  wave  speeds  are  possible  only  if  there  is 
a  sufficiently  large  temperature  gradient  present  in  the  medium  through 
which  the  wave  is  to  propagate.  Qualitative  support  for  this  phenomenon 
based  on  predictions  from  microscopic  theories  of  heat  transport  are 
cited  from  the  literature. 

Finally,  motivated  partially  by  the  results  of  section  5,  we  derive 
in  section  6  a  linear  heat  conduction  equation  governing  infinitesimal 
time-dependent  temperature  variations  superposed  on  a  finite  nonuniform 

a 

equilibrium  temperature  change  from  the  reference  temperature.  The  non¬ 
uniformity  of  the  equilibrium  temperature  change  is  essential  here;  the 
resulting  linear  equation  can  give  rise  to  finite  wave  speeds  and  appears 
to  explain  the  "second  sound"  phenomenon.  The  predicted  thermal  waves 
are  such  that  they  may  travel  in  only  one  direction  relative  to  the 
temperature  gradient  -  either  with  or  against  -  depending  on  the  sign  of 
a  certain  material  constant.  This  is  similar  to  the  relation  between  the 
relative  directions  of  the  heat  flux  and  temperature  gradient  vectors. 


6. 


2 .  Preliminaries .  Rate-independent  response 


Let  x  be  the  (three  dimensional)  position  vector  relative  to  a 
fixed  rectangular  Cartesian  coordinate  system  and  let  t  denote  time. 

Since  we  shall  be  concerned  only  with  rigid  stationary  conductors,  no 
distinction  between  material  particles  and  their  positions  in  space  is 
necessary.  Let  <,  T),  0  and  q  denote  the  internal  energy  per  unit 
mass,  the  entropy  per  unit  mass,  the  temperature  (assumed  positive),  and 
the  heat  flux  vector,  all  functions  of  x,  t  in  a  prescribed  domain. 

These  quantities,  together  with  the  heat  supply  per  unit  mass  r,  must 
satisfy  for  each  x,  t  the  energy  equation-*" 

-  div  q  +  pr  -  pc  ,  (2.1) 

and  the  entropy  production  inequality  which,  for  later  convenience,  we 
write  as 


v  ■  Vi  +  v  £  0  , 

Tloc.  Tcon.  * 


^loc.-O^e^a-of  -  "con.  ■  -  7  a.  ■  &  > 


(2.2) 


where  vloc  and  YCQn  represent  the  local  entropy  production  and  the 
production  of  entropy  due  to  conduction,  the  abbreviation  g  ■  grad  0 
stands  for  the  temperature  gradient,  a  superposed  dot  denotes  partial 


+These  local  forms  can  be  deduced  from  the  corresponding  statements  in 
integral  form.  See,  for  example  Truesdell  and  Noll  [lO] ;  our  q,  r 
correspond  to  their  -h,  q  . 


7 


differentiation  with  respect  to  t  ,  and  p  is  the  mass  density,  a 
function  of  x  only  for  rigid  stationary  bodies.  Combining  (2.1)  and 
(2.2),  we  obtain  the  inequality 


-  pe  +  pQTl  -  |  q  •  g  5  0  .  (2.3) 

In  terms  of  the  free  energy  function  jf  defined  by 

♦  -  «  -  1)0  ,  (2.4) 

the  energy  equation  (2.1)  and  the  inequality  (2.3)  can  be  expressed  in 
the  alternative  forms 

-  div  q  +  pr  ■  p(i  +  1)0  +  Tie)  ,  (2.5) 


and 


-  pi  -  pH0  -  |  •  g  £  0  .  (2.6) 

Taking  as  independent  variables  9  and  g  we  consider  in  this 
section  a  material  whose  response  is  characterized  by  rate-independent 
constitutive  assumptions  of  the  form++ 

t  -  t(e,g)  ,  T)  -  T)(e,g)  ,  q  -  q(e,g)  .  (2.7) 

Lbie  same  symbols  can  be  used  here  for  a  function  and  its  value  without 
confusion. 


8. 


.  -f**vr  '$  *  !V; .  •  ■'  ■ 


The  introduction  of  (2.7)  into  (2.6)  results  in 


-p(§£  +  *n)e  -  p(||) 


g 


g  5  o 


(2.8) 


With  0,  g,  0  fixeu,  all  terms  in  (2.8)  are  determined  except  g  , 
which  may  assume  arbitrary  values.  Then,  in  order  to  satisfy  (2.8),  we 
must  have 


♦  -  t(e) 


(2-9) 


and  (2.8)  reduces  to 


-p(ff  +  H)0  -  •  P  0  .  (2.10) 

Using  (2.10)  and  repeating  the  above  argument  with  0,  g  fixed  and  0 
arbitrary  we  find 

n  *  -  ff  *  •  -  ?  a  •  & 8  0  •  <2-u> 

Thus  the  entropy  is  determined  by  the  free  energy,  a  function  of 
temperature  only,  and  according  to  (2.11)2  heat  cannot  flow  in  a 
direction  of  increasing  temperature. 

Let  the  symmetry  group  of  the  material  be  defined  by  the  set 
of  all  time-dependent  orthogonal  second  order  tensors  A  .  Then,  under 

A/ 

a  change  of  reference  frame  characterized  by  A  ,  the  scalars  e,  Tl,  0, 
if  are  unaffected  but  the  vectors  g  transform  according  to 


9. 


q-Aq  ,  g  -  A  g 


(2.12) 


The  constitutive  equations  (2.7)  are  invariant  under  such  a  change  of 
frame  if  and  only  if  for  each  A  and  all  values  of  q,  g  , 


A  q(6,g)  =  q(0,A  g) 


(2.13) 


If,  in  particular,  the  material  has  a  center  of  symmetry  the  identity 
(2.13)  must  be  satisfied  also  for  A  ®  -1  ,  where  1  is  the  unit  tensor, 
so  that 


-q(e,g)  «  q(e,-g)  . 


(2.14) 


It  follows  that  if  q  is  a  continuous  function  of  g  at  g  *  0  ,  then 


3(e,o)  -  o 


(2.15) 


* 

i.e.,  there  can  be  no  heat  flow  in  the  absence  of  a  temperature  gradient. 

From  the  energy  equation  (2.5)  with  (2.7),  (2.9)  and  (2.1l)  there 
follows  in  component  form 


• 

+  ^ae“^e,k  “  p6^^e  =  pr 


(2.16) 


where  q,  ,  g.(s  0  ),  are  the  Cartesian  components  of  q,  g,  a  comma 

K  X  >  X  ***  *** 


Similar  arguments  have  appeared  many  times  in  the  literature.  See, 
for  example  [2], 
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preceding  an  index  denotes  partial  differentiation  with  respect  to 
rectangular  Cartesian  coordinates  (k=  1,2,3)  >  and  the  usual  summa¬ 
tion  is  implied  by  repeated  indices.  When  q,  *  in  (2.7)-,  (2.9)0  are 

Arf  '  3  C. 

given  explicit  representations  in  terms  of  0,  g  ,  (2.16)  provides  a 
quasi-linear  partial  differential  equation  to  be  satisfied  by  0  .  From 

t 

(2.16),  together  with  (2.15),  it  follows  that  when  the  heat  supply  r 
vanishes  in  rigid  conductors  having  the  rate-independent  type  thermal 
response  (2.7)  there  can  be  no  time  dependence  in  the  temperature  field 
without  a  corresponding  spatial  dependence^  i.e.,  there  can  be  no  change 
in  temperature  except  that  resulting  from  the  flow  of  heat.  Correspondingly, 
in  view  of  (2.1l)2  and  (2.2)^,  the  entropy  production  is  due  to  conduction 
alone  and  there  is  no  local  entropy  production. 

The  residual  energy  equation  (2.16)  is  second  order  in  its  spatial 
derivatives  and  only  first  order  in  its  time  derivative  and  therefore 
clearly  cannot  predict  thermal  waves  (propagating  second  order  discon¬ 
tinuities)  with  real  finite  wave  speeds.  This  will  be  elaborated  upon 
in  sections  4-6. 

We  close  this  section  by  noting  that  the  classical  linear  Fourier 
heat  conduction  equation  results  from  the  linearization  of  (2.16)  when 
(2.7)3  takes  the  special  form  represented  by  the  Fourier  law 

q(e,g)  =  -  K(0)g  (2.17) 

and  ijf(0)  in  (2.9)2  is  assumed  to  be  polynomial. 

_  — 

Assuming,  of  course,  that  does  not  vanish. 

ae 
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We  consider  here  a  generalization  of  the  constitutive  assumptions 
(2.7)  by  including  the  first  time  derivative  of  temperature,  as  well  as 
the  temperature  and  its  gradient,  as  independent  variables .  Thus, 
instead  of  (2.7),  we  now  assume 

*  -  if(e,g,e)  ,  Tj  ■  T|(0,g,0)  ,  q.=  q(e,g,e)  ,  (3.1) 

so  that  the  inequality  (2.6)  becomes 

-p(||+  H)e-p(|*)  •  g-p(7|)e-|  q-gS°  •  (3.2) 

With  0,  g,  8  fixed  it  follows  from  (3*l)  that  all  terms  in  (3.2)  are 

•• 

determined  except  0  ,  which  may  assume  arbitrary  values.  In  order  for 

(3.2)  to  be  satisifed  for  all  0  ,  its  coefficient  (^)  must  vanish,  i.e., 

90 

4  must  be  independent  of  0  .  Similarly  4  cannot  depend  on  g  ,  so 

T  T  f** 

that  (3.1)-^  and  (3.2)  reduce  to 

(i  *  t(e)  (3.3) 

and 

-p(f|+D)e  -  is  •  S5  0  ’ 

respectively. 
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Next  we  decompose  Tj  into  two  parts 


i(e) 


* 

defined  by 


Ti(e,g,0)  -  ,n(o)(0)  +  Ti(e)(e,g,0)  ,  t/o)(0)  =  tK0,o,o)  .  (3.5) 

m  ^  M 


*  * 

Assuming  that  T|  is  continuous  in  0  at  0=0  ,  we  have  in  the  limit" 
as  a  -»  0 


Tl(e,g»a0)  =  11(6,8,0)  +  o(l) 


(3.6) 


so  that 


T^e)(ft,O5a0)  =  o(l) 


(3.7) 


where  a  is  a  real  number.  Now  in  (3*4)  put  g  *  0  ,  replace  0  by 
o/Q  and  use  (3.5)  ,  (3-7)  to  obtain 


-p(||  +  Tl^°^)of0  +  o(a)  5  0 


(3.8) 


where  lim 
cr-0 


0 


It  then  follows  from  (3*8)  that 


(°)  _  .  M 
a© 


(3.9) 


It  shouLd  be  noted  that  the  decomposition  of  Tl  in  (3*5)  does  not  imply 

( o) 

that  T]  is  independent  of  t.  Here  0,  g,  0  are  considered  as  independent 

/  Q  \  ^  s 

variables.  Alternatively  T y  '(0)  is  the  limit  of  T)(0,ae,60)  as  a,0  -*  0  . 

A  function  of  a  defined  in  a  neighborhood  of  is  o(l)  or  0(l)  as 

a  -*  aQ  according  as  it  is  zero  or  bounded  in  this  limit. 
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and  the  inequality  (3*4)  assumes  the  form 


-  pH(e)e  -  *  S  £  0  •  (3.:.o) 

In  view  of  (3.1),  (3*3),  (3*5)  and  (3-9)  the  energy  equation  (2.5) 
becomes  (in  component  form) 


which  is  the  counterpart  of  (2.16)  under  the  assumptions  (3.1)* 

Many  of  the  conclusions  drawn  in  section  2  are  altered  for  the  rate- 
dependent  response  under  consideration.  In  particular,  even  though  i|r 
still  reduces  to  a  function  of  0  only,  we  cannot  conclude  the  same 
about  T1  which  remains  dependent  on  g  and  0  .  Furthermore,  it  is 
clear  from  (2.2)  and  (3.10)  that  the  inclusion  of  0  as  an  independent 
variable  in  constitutive  assumptions  (3.1)  has  led  to  a  local  entropy 
production.  In  addition,  the  total  entropy  is  no  longer  expressible  in 
terms  of  the  free  energy  as  in  the  case  not  only  of  the  rate-independent 

+It  is  worth  noting  that  if  in  place  of  (3.1)  we  had  assumed  ^(0,g,0,g), 

etc.,  that  is,  if  the  rate  of  the  temperature  gradient  had  been  in-  ~ 
eluded  as  an  independent  variable  also,  then  (3-3)  would  be  replaced  t>y 
This  should  be  contrasted  with  the  results  obtained  by 
Coleman  and  Mizel  [ll3,  where  is  it  shown  that  t  cannot  depend  on 
temperature  gradients  when  a  dependence  on  rates  is  excluded.  See 
also  in  this  regard  Eringen  C 12 J  where  the  rate  of  temperature  and 
rates  of  temperature  gradients  are  included  as  independent  variables 
in  the  context  of  a  general  thermo-mechanical  investigation. 
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type  response  but  also  of  the  more  general  response  characterized  by  a 
functional  over  the  time  history  of  temperature  as  investigated  by 
Coleman  and  Gurtin  [43. 

The  requirements  that  the  constitutive  equations  (3.3)  and 
(3.1)9  ,  satisfy  the  appropriate  transformation  relations  under  a  change 
of  reference  frame  (imposed  by  the  symmetry  group  of  the  material),  that 
the  heat  flux  function  q  be  continuous  in  g  at  g  =  0  and  that  the 

/NV  /V  A# 

material  possess  a  center  of  symmetry  imply  (with  an  argument  parallel 
to  that  used  to  arrive  at  (2.15))  that  q  must  be  an  odd,  and  T|  an 
even,  function  of  g  with 


q(e,O,0)  =  0  .  (3.12) 

Therefore,  as  in  the  rate -independent  case,  there  can  be  no  heat  flux 
corresponding  to  zero  temperature  gradient.  We  also  observe  the 
inequalities 


q(e,g,o)  •  g  $  o  ,  n(e)(0,o,e)e  *  o  ,  (3.13) 

/■w  *** 

which  follow  from  (3.10),  but  we  note  that  q  *  g  5  0  cannot  in  general 
be  concluded. 
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4.  Propagation  of  thermal  waves 


We  investigate  here  the  possibilities  of  the  existence  of  propagating 
second  order  discontinuities,  i.e.,  discontinuities  in  the  second  deriva¬ 
tives,  in  the  temperature  field.  For  convenience,  however,  we  confine 
our  analysis  to  one  space  dimension  x  .  Writing  q,  g  for  the  x- 
components  of  q,  g  ,  we  obtain  from  (3*11)  its  one-dimensional  counter- 
part 


<H)e  „  +  (*?  *  P0  ^)e  „  ♦  p(e  ^)e 


^g,v,xx 


ae 


ag  /v,x 


ae 


+  <>,x  *  ^  +  0  ^  -  6  ’  0r 

00 


(4.1) 


We  assume  that  0,  g,  9  and  r  are  continuous  functions  of  x,  t 
and  denote  by  [  ]  the  value  at  the  propagating  discontinuity  of  the  jump 
of  whatever  quantity  it  brackets.  If  we  allow  discontinuities  in  second 
order  derivatives  of  0  and  assume  that  the  first  derivatives  are 
continuous,  there  follows  from  (4.1) 


‘S'-lr 


)£e  J  +  p(e 

,X  30 


=  0  . 


Using  standard  jump  conditions 


we  have 


(4.2) 


* 


See  Thomas  C 133  or  Truesdell  and  Toupin  [l4;  Sec.  l8l). 
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(4.3) 


uCe  ] 


XX 


[§)  =  iftej 


where  U  is  the  speed  of  propagation  of  the  discontinuity.  Substitution 
of  (4.3)  into  (4.2)  gives 


2a  +  „e  stel  t  [(aa  +  pe  2I&2  -  4p6  |a 

ae  as  ae  36  ag  ae 


(e) 


M. 


u  = 


2p0 


3g 

TV 


(4.4) 


an 


ae 


i  ^n(e)  *n(e  ^ 

provided  [0  -i  and  - 1  ■■  are  not  zero.  In  the  event  — 1 -  does 

,xx  ae  ae 

vanish, (4.2)  and  (4.3)  yield 


U  = 


& 


&  +  peM 


TV 


(4.5) 


ae 


ag 


The  terms  on  the  right  hand  side  of  (4.4),  (4.5)  are  evaluated  at  the 
propagating  discontinuity  so  that  U  can  be  considered  alternatively 
as  a  function -of  x  or  t  only. 

The  following  observations  are  apparent  from  (4.4)  and  (4.5): 

(i)  Equation  (4.4)  shows  that  second  order  discontinuities  in  temperature, 
i.e.,  thermal  waves,  may  possibly  occur  with  two  distinct  real  wave 


speeds.  This  depends,  of  course,  on  the  relative  magnitudes  and  signs 
of  the  terms  in  (4.4)  and  will  be  investigated  further  for  particular 


An  alternate  derivation  of  (4.4)  results  from  the  characteristic 
relation  appropriate  to  the  quasi-linear  partial  differential  equation 
(4.1).  It  is  thereby  seen  that  discontinuities  of  the  type  under 
consideration  can  occur  only  along  a  characteristic.  See  Courant 
and  Hilbert  C 15  J . 
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(e) 

representations  of  Ip  7 ,  if  and  q  in  the  following  sections. 

(ii)  It  follows  from  (4.5)  that  if  is  independent  of  6  ,  then 

only  one  wave  speed  is  possible.  In  fact  we  see  that  thermal  waves  may 

(g) 

occur  when  Tp  vanishes  altogether.  Thus,  thermal  dissipation  due 
to  a  local  entropy  production,  represented  by  the  first  term  of  the 
inequality  in  (3.10),  may  not  be  necessary  for  thermal  waves,  provided 
q  remains  dependent  on  0  as  well  as  g  .  Equation  (4.5)  also  siiowB 

(g]  • 

that  when  Tp  7  vanishes  and  q  does  not  depend  on  0  ,  as  In  the 
rate- independent  response  discussed  in  section  2,  no  finite  Rperi 
thermal  wave  can  occur. 
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5 .  Linearity  in  temperature  rate 


Recalling  (3.1),,  ,  and  (3-3)  we  introduce  in  this  section  specific 
assumptions  regarding  the  dependence  of  the  response  functions  on  the 
temperature  rate.  Before  doing  so  we  decompose  g  (as  was  done  with 
in  (3.5))  into  two  parts  defined  by 

^(e.g.e)  -  </°^(®.jg)  ♦  q^(s.jg.®)  »  j/°)(®>jg)*;j(®>g>0)  •  (5*1) 

>ien  (LU)  l*i>llea 

i ( i  ,0 , S )  -  i(o,(S.O)  ♦  q<e)(0,O,e)  -  0  .  (5.2) 

4  «  4  M  MS  MS 

*h«  )  an l she a  t»  (5.1),  it  follows  that 

i,#'.  ’’  •  J (oM*.0)  -  £  •  (5.3) 

I  i«  (5. I)  we  have 

^"'(0.0)  -  0  ,  ^e)(«,£,0)  -  0  .  (5.4) 

l#nH  f  ui  I  |»  lentiy  V ,  are  arbitrary  functions  of  0  , 

I  1  -ii  .  nitlniHiua  fut  ‘tion  of  0,  g  which  is  odd  in  g 
i  *i (  4)  Sr  aaaumr  that  the  "extra"  parts  of  entropy  and 

«i  ii  i>  namely  ^  and  ^ .  are  nonlinear  functions  of  0,  g 

a 

an*  a r*  linear  fun  t Iona  of  degree  one  in  0  so  that 

V'Nl.g.t)  •  -  n(0,g)0  ,  *  *(0»g)®  •  (5*5) 
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Assuming  further  that  q^,  q^  are  polynomials  in  g  and  recalling 
(5.4),  we  arrive  at  the  forms 

q^(e,g)  -  -  K(e,g)g  ,  q^(e,g,e)  =  H(e,g)ge  .  (5.6) 

^  M  MM  MW  W  MW  MW  MS 

By  (3.13),  (5*1)  and  (5*2)  the  functions  n(0,g),  K(e,g)  in  (5*5), 

MW  MW  M/ 

(5.6)  must  satisfy 


n(9,0)  *  0  ,  K(0,g)g  •  g  5  0  . 

MW  MW  MW  MW  MS 


(5-7) 


In  view  of  (5-6) ,  the  energy  equation  (3. 11)  now  assumes  the  form 


["  agj  ^  gP  +  agx  g^®]0,kjt  +[Hkj  gj" pe®  agk]V 

-  cne5  + 


-  p(ne  +  ee  +  e  M)e  =»  pr  , 
90  36 


(5.8) 


Which  is  a  quasi-linear  second  order  partial  differential  equation  to  be 
satisfied  by  the  temperature.  The  classification  and  type  of  initial¬ 
boundary  value  problem  appropriate  to  such  an  equation  has  been  exten- 

* 

sivaly  discussed  in  the  literature. 

The  essential  character  of  equation  (5-8)  will  not  be  altered  if  we 
make  some  further  simplifying  assumptions.  In  particular,  if  we  assume 


* 


See  for  example  [  15  ] . 
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that  n,  K  and  H  introduced  in  (5-5)  and  (5*6)  are  functions  of  8 
only,  then  (5-8)  becomes 


(■\j  +  “kj  ^0,kj  +  “kj  gj  ®,k 

ax.  .  a& 

-  pnee  +  +  (“ae 

-  p(n8  +  88  ~  +  0  -^-^)8  =  pr  .  (5*9) 

36  ae 

If,  in  addition,  we  assume  the  material  is  isotr'  oic  with  a  center  of 
symmetry,  i.e., 


\j(0)  *  >  Ve)  • 


(5.10) 


then  (5-9)  becomes 


(-k+he)8>kk  +  h  ^  0jk  -  pn08 

+  ^  ae  gk  +  38  6  gk  ® ,k 

2 

-  p{ne  +  80  |r  +  0  ^-|)9  *  pr 

36  a  9^ 

and,  in  view  of  (5 .7)j 


(5.11) 


k  5  0  .  (5.12) 

We  now  use  the  one-dimensional  form  of  (5>ll)  to  examine  more 
closely  the  propagation  of  thermal  waves.  Following  the  procedure 
that  led  to  (4.4)  we  then  obtain  for  the  speed  of  propagation  U  of 
second  order  discontinuities  the  expression 
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-  (h0  v  t  [(he  J2  -  4pne(k-  he)]1/2}/2pne  ,  if  n  i  0  , 

>x  )x 


U  -  <  (he  -  k)/he  v  ,  if  n  -  0  ,  he  ^  0  (5.13) 

I  ’  ’ 

'  oo,ifn*0,h6  *0  (classical  case) 

Equation  ( 5 .13) ^  reveals  an  interesting  phenomenon.  No  real  wave  speed 
can  result  if  both  e  and  6  are  zero  "in  front  of  the  wave,"  that 

ix 

is,  if  the  conductor  into  which  the  wave  is  to  propagate  is  in  a  uniform 
equilibrium  (time- independent)  thermal  state.  On  the  other  hand,  real 
wave  speeds  can  clearly  occur  provided  the  temperature  gradient  in  the 
conductor  is  large  enough.  This  result  appears  to  be  in  harmony  with 
the  proposals  made  by  Chester  [5]  and  others  on  a  microscopic  basis  that 
thermal  waves  should  be  detectable  in  many  solids  but  only  above  a 
certain  critical  state  of  thermal  agitation. 
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6 .  Infinitesimal  time-dependent  temperature  variations  superposed  on  a 
finite  equilibrium  temperature  change  from  reference 

It  was  found  in  the  previous  section  that  a  thermal  wave  cannot 
propagate  into  a  rigid  conductor  with  a  rate-dependent  thermal  response 
unless  the  magnitude  of  the  temperature  gradient  is  sufficiently  lai/e . 
This  fact  suggests  the  possibility  of  predicting  thermal  waves  charac¬ 
terized  by  a  linear  differential  equation  governing  infinitesimal  time- 
dependent  temperature  variations,  provided  we  superpose  these  time- 
dependent  variations  on  a  finite  time-independent  temperature  change 
from  the  constant  reference  temperature.  Here  for  simplicity  we 
confine  our  attention  to  an  isotropic  material  and  use  (5.11)  . 

Let 


e(x,t)  =  e(x)  +  €0'(x,t)  ,  r(x,t)  =  r(x)  +  er'(x,t)  ,  (6.1) 

^  A/  ^  M  ^  M 

i.e.,  we  consider  temperature  variations  from  equilibrium  that  are  of 
0(e)  in  the  limit  as  +  r.  —  0.  Also  let 

e(x)  =  e  +  e(x)  ,  (6.2) 

A 

so  that  0(x),  which  is  not  necessarily  uniform,  represents  the  finite  time- 

independent  variation  from  the  reference  temperature  @q.  The  temperature 
A  . 

variations  0  and  0  may  be  positive  or  negative  within  the  limitations 

+In  this  section  e  is  a  small  dimensionless  parameter.  It  is  implicit 
in  (6.1)  that  the  time  and  spatial  derivatives  of  the  temperature 
variation  9  -  is  also  of  0(e)  . 
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imposed  by  (6.1)^,  (6.2)  and  the  assumption  0  >  0  . 

Corresponding  to  (6.1)  we  introduce  the  following  decomposition  of 
entropy  T|  through  its  parts  1^°^ ,  T^e),  namely 


H(o)  -  n(o)(e)  -  ii(o)(x)  *  en'<o)(x,t)  , 

T1(e)  =  Tl(e)(e,g,0)  =  Tl(e)(x)  +  eT1/(e)(x,t) 

t**  r+s 


(6.3) 


From  (5*5)^,  (6.1)  and  (6.3)g  we  find 

^j(e)  =  0  5  ^/(e)  =  .  }  (6.4) 

so  that  the  entropy  variation  from  its  equilibrium  value  T^°^(x)  is  also 

'(V 

Of  o(e)  . 

Next  we  denote  by  \lr(x)  the  value  of  the  free  energy  at  the 
equilibrium  state,  i.e., 


♦  (x)  =  f[e(x)] 

/sw  rv/ 


(6.5) 


and  let  4f'(x,t)  stand  for  the  time-dependent  variation  of  ^  from 

if  so  that 


%'(x,t)  =  ty(x,t)  ’  t(x)  •  (6-6) 

Consistent  with  time -dependent  variations  in  temperature  and  entropy 

2 

that  are  of  0(e),  ^'(x,t)  must  be  of  0(e  )  and  it  is  sufficient  to 
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assume 


++ 

^'(x,t)  =  -  h  c[ee'(x,t)]2 

C.  #"W 


(6.7) 


but  we  make  no  specific  assumption  regarding  the  dependence  of  ty  on 
9  in  (6.5).  Noting  from  (6.1)^  that 


_3_ 

30 


_a  1  a 

30  6  30' 


3 


we  obtain  from  ( 3 -9) s  (6.5) »  (6.6)  and  (6.7) 


(6.8) 


T1(°)  =  .  M+  c(c0')  ,  (6.9) 

30 

and  hence  by  (6.3)  there  results 


~(°)  M  „,(o) 

7)'  - - £  ,  V  =  c0 

30 

Also,  we  have 

2  2~ 

Ui  =  UL  .  c 
2  ~2 

30  30 


(6.10) 


(6.11) 


Making  use  of  (6.11)  and  (6.1)  in  (5. 11),  for  an  isotropic  rigid 
conductor  we  obtain 

++We  could  include  in  (6.7)  a  linear  term  c..e9'  but  it  enters  the 

~(o) 

development  only  through  its  coefficient  c^  as  a  constant  in  T|  '  and 
hence  can  be  absorbed  in  the  value  of  T|  at  the  reference  temperature  0q  . 
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(6.12) 


*1  e'u,  6>kk  «'♦*  9>kk  «'  +  h  ejk  e;k 


r**/  • 


M  m  w  ^  /v  m  <vO  ^  <vO  • 

-p  n  e  e;  -  kx  e^k  e'k  -  k2  e  k  e  +  hx  e  k  e' 


c)Q'  -  pr'}  +  0(e2)  =  0  , 

where  k,  k^,  kg,  h,  and  n  are  defined  by  the  Taylor  expansions 


k(9) 

*  k+  kx  c9 

+  I^Ue')2  h 

-  0(€3) 

,  n(9)  =  n  +  0(e) 

m 

■  h  4-  h,  *9 ' 

+  0(€2)  ;  k 

-  k(9) 

k  =  —  1 
’  1  d9  ~ 

0  1  9=9 

** 

*2  “ 

d2k 

~2 

y 

n  »  n(9)  , 

h  =  h(9) 

~  dh  . 

*  hl  "  d9  L  ~ 

d9 

/V 

e-9 

|  9=9 

assumed  to  be  valid. 

Equation  (6.12)  is  satisfied  to  0(e)  if  tbj  finite  time-indepen- 

f+J 

dent  variation  9(x)  satisfies 


k 


Pr  » 


2 

and  to  0(e  )  if  in  addition  to  (6.l4)  the  infinitesimal  time  dependent 
variation  9,(x,t)  satisfies 


-k  e'kk  +  h  9>k  9'k  -  p  n  0  9' 


kl  0 ,k  0*k  +  0,kk 


/-wp 

+  hn  9 


»k 


4 


P*<«  -  4i>lr 

ae 


(6.13) 


(6.14) 


(6.15) 
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Since  k,  k^,  k^,  n,  h  and  depend  only  on  0  ,  the  coefficients  in 
(6.15)  are  determined  once  the  solution  of  (6.14)  is  known. 

Equations  (6.14)  and  (6.16)  determine  the  infinitesimal  time- 
dependent  temperature  variation  superposed  on  a  finite  equilibrium 
change  from  the  reference  temperature  when  the  conducting  medium  is 
characterized  by  the  rate-dependent  constitutive  relations  of  sections 
3  and  5.  The  corresponding  results  for  the  rate-independent  response 
discussed  in  section  2,  for  an  isotropic  medium,  can  be  obtained  if  we 
put  n  =  h^  =  h^  =  0  .  Then,  (6.15)  becomes 


-k  0' 


,kk 


kl  9,k 


9,,.+  p6(c-  ^4)0' 
ae 


,k 


<kl  0,kktk2  6,k)e  *  l“r 


where  0  continues  to  satisfy  (6.l4).  If,  in  particular,  0 

coincides  with  0  ,  then  0'  is  measured  from  the  uniform  reference 

o 

temperature  Qq  ,  ^  will  have  a  value  ,  r  must  vanish  according 
to  (6.l4),  and  (6.l6)  reduces  to  the  classical  form 


-k  e'kk  *  f  t  «a  i'  ■  Pr’  , 

/S/ 

where  k  ,  c  0Q  correspond  to  the  usual  thermal  conductivity  and 
specific  heat. 

The  generalization  of  the  heat  conduction  equation  represented  by 


(6.16) 


(6.17) 
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(6.15)  is  &  linear  second  order  partial  differential  equation  with  coef' 


~  * 
ficients  depending  only  on  x  through  0  and  is  in  standard  form. 

Its  type  (elliptic,  parabolic,  or  hyperbolic)  and  hence  the  associated 

initial -boundary  value  problem  that  is  well  posed  depends  on  the  signs  of 

its  coefficients  when  reduced  to  canonical  quadratic  form.  In  particular, 

^  ^  /v 

when  h,  n,  h.  vanish  it  reduces  to  the  classical  parabolic  equation.  When 
these  coefficients  do  not  vanish  but  0  .  vanishes,  i.e.,  when  the 
equilibrium  temperature  field  is  uniform,  then  the  equation  is  elliptic 
since  k,  p  n  0  are  both  non-negative.  No  finite  speed  waves  are 
possible  in  either  of  these  circumstances . 

The  complete  classification  of  (6.15)  for  the  one  space  dimensional 
case  is 


-  Up  k  n  0 


<  0  ...  elliptic 
*  0  . . .  parabolic 
>  0  . . .  hyperbolic 


(6.18) 


The  associated  speeds  of  propagation  U  of  second  order  discontinuities 
are 

M  SS*  O  ^  ^.1  /O  A/  A/ 

1{-  h  0  t  [(h  0  J  -Up  k  n  0]X/H/2P  n  0  ,  n/0 

,x  ,x 

-  k/h  0  ,  n«0  ,  h?x/0  (6.19) 

Arf  A/  A# 

«  ,  n-0  ,  he  *  0 


For  the  classification  of  this  equation  and  a  reduction  to  its 
canonical  form  see  [15>  p.  l8l] . 
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We  see  from  (6 .19)  that  the  wave  speed  U  must  have  the  opposite 


algebraic  sign  from  h  0  .  Hence,  if  h  is  positive,  the  present 

>x 

theory  predicts  that  thermal  waves  can  propagate  only  in  the  direction 
from  hot  to  cold,  the  same  direction  as  the  part  of  the  heat 

flux.  The  entropy  production  inequality  places  no  restriction  on  the 
sign  of  h  .  Fran  (5*6),  (5.10)  we  see  that  the  "extra"  part  of  the 

(e) 

heat  flux  q'  '  changes  its  direction  with  respect  to  the  temperature 
gradient  as  the  sign  of  0  changes. 
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